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Abstract 

Our first-principles computations show that the ground state of 
PbTiOa under hydrostatic pressure transforms discontinuously from 
PAmm to R3c at 9 GPa. Spontaneous polarization decreases with in- 
creasing pressure so that the R3c phase transforms to the centrosym- 
metric R3c phase at around 30 GPa. The first-order phase transition 
between tetragonal and rhombohedral phase is exceptional since there 
is no evidence for a bridging phase. The essential feature of the R3c 
and R3c phases is that they allow the oxygen octahedron to increase 
its volume Vb at the expense of cuboctahedral volume Va around a Pb 
ion. This is further supported by the fact that neither the R3m nor 
Cm phase, which keep the Va/Vb ratio constant, is a ground state 
within the pressure range between and 40 GPa. Thus tetragonal 
strain is dominant up to 9 GPa, whereas at higher pressures efficient 
compression through oxygen octahedra tilting plays the central role 
for PbTiOs. Previously predicted pressure induced colossal enhance- 
ment of piezoelectricity in PbTi03 corresponds to unstable Cm and 
R3m phases. This suggests that the phase instability, in contrast 
to the polarization rotation, is responsible for the large piezoelectric 
properties observed in systems like Pb(Zr,Ti)03 in the vicinity of the 
morphotropic phase boundary. 

In ferroelectric crystals electric polarization occurs spontaneously. The 
spontaneous polarization occurs when the energy gained by the interaction 
of the dipoles induced by the local field is larger than the energy required 
to produce the dipoles pQ. This viewpoint was adopted once early attempts 
for explaining ferroelectricity were developed: the ferroelectric state is sta- 
ble when the attractive dipolar energy exceeds the distortion energy. It soon 
turned out that the model in its simplest form was insufficient. For BaTiC>3 a 
phenomenological model, based on Slater's 'rattling' Ti ion and the fact that 
the environment of the oxygen ions is noncubic, was developed to explain 
ferroelectricity. A rattling Ti ion creates large ionic polarizability (implying 
small distortion energy), whereas the noncubic environment of oxygen ions 
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increases the local field contribution from the titanium. A large local field 
implies large dipolar attraction energy dominating over the distortion en- 
ergy. Quantitative treatment of ferroelectricity in BaTi0 3 and PbTi03 was 
provided by first principles computations which, instead of a rattling tita- 
nium ion, considered the hybridization between Ti 3d and oxygen 2p states 
[2]. These computations were able to explain the difference between the zero 
pressure ground states of BaTi0 3 (rhombohedral) and PbTi0 3 (tetragonal) 
by considering Ba essentially as a closed shell ion and taking the hybridization 
of Pb 6s and oxygen 2p states into account. The ground state is very sensitive 
against primitive cell volume and one may expect significant changes in the 
structure of PbTi0 3 as a function of hydrostatic pressure. Recent develop- 
ment on ferroelectrics are due to the density functional theory (DFT) which 
provides an accurate value for the total energy except for the errors due to 
the approximate forms used to describe the exchange-correlation functional. 

We used the ABINIT DFT code to compute the total energies at different 
pressures [3] . The computations were carried out within the local-density ap- 
proximation (LDA) approximation and a plane wave basis. Norm-conserving 
pseudo-potentials (scalar-relativistic for Pb and non-relativistic for Ti and O 
[I]) were generated using the OPIUM pseudopotential generator package [5]. 
The valence states included 5d6s6p (Pb), 3s3p3dAsAp (Ti) and 2s2p3d (O) 
states. The energy cutoff was 75 hartree and the fc-point meshes were 6x6x6, 
except for the R3c and R3c symmetries, for which it was 4x4x4 (due to 
the larger primitive cell), and IAcm symmetry, for which a mesh 6x6x4 
was used. Denser fc-point meshes were used for computing spontaneous po- 
larization through Berry's phase method. Structural optimization was car- 
ried out till force acting on each atom was below 5 x 10~ 6 eV/A. Phonon 
frequencies and the corresponding atomic displacements (eigenvectors) were 
computed through response function method implemented in ABINIT. To 
tackle the problem of structural optimization we selected various potential 
crystal symmetries for which the enthalpies were computed as a function of 
pressure. Thus, due to the polarization rotation theory [6] the case where 
the Cm symmetry bridges the tetragonal PAmm and rhombohedral R3m 
symmetries is a particularly attractive possibility to be examined. The Cm 
phase is a common subgroup of the PAmm and R3m according to the group- 
subgroup chains PAmm — > Cmm2 — > Cm and R3m — > Cm. In such a case a 
large electromechanical coupling coefficient, assigned to the Cm phase, was 
predicted [TJ. This was based on the idea that the Cm phase would allow 
energetically favourable continuous rotation of polarization between the [001] 
[PAmm phase) and [111] (R3m phase) directions. At higher pressures the 
R3m phase was proposed to transform to a cubic Pm3m phase. By sym- 
metry, these phases do not possess octahedral tiltings: the role of the Cm 
phase is to be an energetically favourable path between the PAmm and R3m 
phases. An alternative path is a commensurate transition involving octahe- 
dral tiltings, which are frequently observed in oxide perovskites. This was 
suggested to occur also in PbTi0 3 under high pressure according to the phase 



transition sequence P4mm — > 1 4cm — > 1 4/ mem — > 1 4cm [8]. The proposed 
octahedral tilt system is the same as is found in SrTi03. However, octahe- 
dral tilting is also found in rhombohedral R3c and R3c phases. The R3c 
space group is a generic rhombohedral symmetry in the sense that it allows 
both octahedral tilting and cation displacements. Thus, it is important to 
consider both R3m symmetry (only cation displacements are allowed) and 
its subgroup R3c. The latter symmetry corresponds to the case where the 
two oxygen octahedra of the primitive cell are tilted about the pseudocubic 
[111] direction. The study of phonon instabilities provides invaluable infor- 
mation: particularly the phonon frequencies at the high symmetry points of 
the Brillouin zone and the acoustical branch provide very useful information 
related to the possible phase transitions. The first mentioned case provides 
information about commensurate phase transition whereas the incommen- 
surate phase is signalled by a dip in the lowest frequency acoustical phonon 
branch. An incommensurate phase often bridges phases which do not possess 
group-subgroup relationship. 

Figure [T] shows the enthalpies for P4mm, Cmm2, Cm, 1 4cm, R3m and 
R3c phases with respect to cubic Pm3m phase. As far as P4mm and R3m 
phases are concerned, it bears a close resemblance with the enthalpy curves 
given in ref. [TJ. The difference between the enthalpies is very small except 
for the rhombohedral symmetries. This is due to the fact that also in the 
case of the Cm and Cmml symmetries the structural parameters turned out 
to correspond to the P4mm though orthorhombic and monoclinic distor- 
tions were allowed. Neither orthorhombic nor monoclinic was energetically 
favourable as Figure [2] shows. Figure [1] reveals that the minimum enthalpy 
path corresponds to the phase transition sequence P4mm — > R3c — > R3c. 
It is worth mentioning that, within the studied pressure range from to 40 
GPa, neither the Cm nor R3m phase was a ground state although with in- 
creasing pressure the Cm phase became more favourable than P4mm phase 
and the R3m phase was more favourable than the Cm phase at even higher 
pressures. For the electromechanical properties it is crucial to note that 
the phase transition between the P4mm and R3c phases must be of first 
order and thus does not involve continuous polarization rotation. To put 
these observations on a firm ground phonon frequencies were computed at 
the high symmetry points of the Brillouin zone as a function of pressure. 
At zero pressure no phonon instabilities in the P4mm phase were observed, 
consistently with experiments. However, once the pressure is increased the 
Brillouin zone corner (A point, q = (---)) B\ mode becomes unstable at 10 
GPa, Figure EJ By further increasing the pressure also the E symmetry mode 
becomes unstable at 13 GPa. No other instabilities or anomalies in acoustical 
branches were observed. Though the frequencies of the Brillouin zone centre 
Al(ITO) and .E(ITO) modes were decreasing with increasing pressure they 
frequencies remained finite at the whole pressure range. Figure H] a shows 
the oxygen displacements in B\ and E symmetry and the structural distor- 
tion following the condensation of both modes. This bears a close analogy 



to the Pm3m — * R3c phase transition observed in LaA103 [9]. The fact 
that the B>i mode becomes unstable first suggests that the phase transition 
may occur in two stages. The oxygen octahedral rotation corresponding to 
the condensation of the B\ mode is depicted in Figure H] b and is similar 
to the case observed in SrTiC>3 [9j. The second stage is shown by black ar- 
rows in Figure H] a. However, the IAcm symmetry did not correspond to the 
minimum enthalpy, Figure [TJ This implies that either the stability range of 
IAcm phase is very narrow or the phase transition between PAmm and R3c 
phase takes place without an intermediate phase. In the case of a first-order 
phase transition the frequencies of the unstable modes do not reach zero at 
the transition point. It is remarkable that the information provided by the 
phonon instabilities, most notably by the phonon eigenvectors, is perfectly 
consistent with the analysis based on enthalpies. 

The Va/Vb ratio decreases with increasing hydrostatic pressure in rhom- 
bohedral R3c and R3c phases. The central role of the tilt angle of the octa- 
hedron about [111] direction (u) is neatly given by the expression Va/Vb ~ 
6 cos 2 uj — 1 [TU] . This is related to the increasing octahedral tilting and is 
clearly an important mechanism for minimizing the enthalpy. This is seen by 
comparing Figures [1] and [5] a and noting that R3m phase does not permit oc- 
tahedral tiltings. Figure [5] b shows that spontaneous polarization decreases 
monotonically with increasing pressure. There is a discontinuous drop in 
spontaneous polarization at 9 GPa so that not only is the polarization di- 
rection changed abruptly but also the magnitude drops significantly. With 
increasing pressure the cation shifts diminish and a transition between the 
R3c and R3c phases occurs at around 30 GPa, see Figure [TJ At 30 GPa the 
cation shifts from the centrosymmetric positions are small, which suggests 
that the transition is continuous. These results mean that under high pres- 
sure PbTi03 does not minimize its enthalpy through polarization rotation 
but by tilting the oxygen octahedra. This has very important consequence 
for the piezoelectric properties: though the di§ constant has a clear maxi- 
mum at the phase transition there is no dramatic enhancement of piezoelec- 
tric constants at the phase transition pressures, as Figure c shows. This 
is in complete contrast to the predictions based on the polarization rota- 
tion theory [6], [7] : in the pressure range between 9 and 15 GPa significantly 
larger piezoelectric coefficients are seen in the case of the R3m and PAmm 
phases, which do not permit polarization rotation. We conclude that it is 
the phase instability which results in the large piezoelectric constants. In 
a classical piezoelectric Pb(Zr,Ti)03 system an analogous situation occurs 
in the vicinity of the morphotropic phase boundary (MPB) where the Cm 
phase coexists with the R3m phase (at room temperature) or R3c phase (at 
low temperature) [UJ . Our computations suggest that it is the phase insta- 
bilities of the R3m and/or Cm phase(s) which is responsible for the large 
electromechanical coupling coefficients observed at MPB. 

Experimentally this zero temperature phase transition sequence PAmm — > 
R3c —>■ R3c can be verified through neutron diffraction experiments. Neutron 



diffraction allows oxygen positions to be precisely located and the superlattice 
reffections allow the identification of octahedral tiltings. In addition, Raman 
spectroscopy provides useful complementary information: all phases have 
very distinctive spectra. Since the transition between R3c and R3c is likely 
to be continuous, Raman scattering data are especially valuable: the Bril- 
louin zone centre normal modes (IR and R stand for the infrared and Raman 
active modes) of R3c and R3c transform as the representations 5Ai(IR, R) © 
A 2 © 10E(IR, R)@ and A 2g © E g (R) © AA lu © hA 2u (IR) © 9E U (IR), respec- 
tively. Thus, the R3c phase has 15 peaks (excluding accidental degeneracies) 
whereas R3c has only one peak. 
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Figure 1: Enthalpies as a function of pressure with respect to the cubic 
Pm3m phase. Our first-principles computations reveal that energetically 
favourable and large crystal compression is allowed by tilting the oxygen 
octahedra. The minimum enthalpy path is PAmm — > R3c — > R3c. 
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Figure 2: Lattice parameters for the P4mm (crosses) and Cm (open circles) 
phases. For an easier comparison, the monoclinic a and b axes are divided 
by y/2. Inset shows the behaviour of the monoclinic j3 angle as a function 
of pressure. It is seen that monoclinic distortion is essentially zero in the 
studied pressure range and thus PbTiOs has PAmm symmetry up to 9 GPa 
at which a transition to the R3c phase takes place. 
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Figure 3: Brillouin zone boundary modes at A point (q 
unstable at 10 GPa [B\ mode) and at 13 GPa (E symmetry mode). Lines 
are guides for the eyes. 



Figure 4: Phase transition between P4mm and R3c phases corresponds to 
the condensation of two modes, a, Oxygen displacements in two unstable 
modes corresponding to the Brillouin zone boundary A point. White and 
black arrows indicate the oxygen displacements in the B\ and E symmetry 
modes, respectively. Grey arrows indicate the superposition of the B\ and E 
symmetry modes resulting in rhombohedral symmetry (the shown octahedra 
is rotated clockwise and the next oxygen octahedra (not shown) along the 
[111] direction is rotated anticlockwise). As a result of the condensation of 
these modes the symmetry discontinuously changes to a rhombohedral R3c 
symmetry, b, Oxygen octahedral tilting in I4cm phase. B\ symmetry mode 
condensation corresponds to the phase transition between PAmm symmetry 
and its subgroup 14cm. The ab basal plane of the 14cm unit cell is indicated 
by dark grey shading. 
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Figure 5: a, The ratio between cuboctahedra (Va) and oxygen octahedra 
volume (Vb) as a function of pressure in R3c (filled marks) and R3c (open 
marks) phases. The corresponding octahedral tilting angle u is also shown, 
b, Spontaneous electric polarization across the proposed phase transition se- 
quence, c, Absolute values of the piezoelectric constants for the stable phases 
(continuous line with filled symbols) and for the unstable phases (P4mm: 
dashed lines, R3m: dashed-dotted line). Squares, diamonds, triangles and 
spheres indicate d 3 i, d 33 , ci 15 and dn constants, respectively. The anoma- 
lously large di$ values of the R3m and P4mm phases at 15 GPa pressure 
are a clear indication of a structural instability Symmetry dictates that for 
these phases polarization rotation is not possible and thus cannot offer an 
explanation for the anomaly. 



